Loschmidt echo in a system of interacting electrons 
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We study the Loschmidt echo for a system of electrons interacting through mean-field Coulomb 
forces. The electron gas is modeled by a self-consistent set of hydrodynamic equations. It is observed 
that the quantum fidelity drops abruptly after a time that is proportional to the logarithm of the 
perturbation amplitude. The fidelity drop is related to the breakdown of the symmetry properties 
of the wave function. 
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Introduction. — In a famous controversy with L. Boltz- 
mann, J. Loschmidt pointed out that, if one reverses the 
velocities of all particles in a physical system, the latter 
would evolve back to its initial state, thus violating the 
second law of thermodynamics. The main objection to 
this argument is that velocity reversal is a very unstable 
operation and tiny errors in the reversal quickly restore 
normal entropy increase. 

More recently, the original idea of Loschmidt was re- 
vived in the context of quantum information theory. In- 
deed, any attempt at coding information using quantum 
bits is prone to failure if a small coupling to an uncontrol- 
lable environment destroys the unitary evolution of the 
wave function (decoherence) . In order to estimate the 
robustness of the system against perturbations from the 
environment, the following procedure has been suggested. 
The system is allowed to evolve under the action of an 
unperturbed Hamiltonian until time T; then it is evolved 
backwards in time until 2T with the original Hamilto- 
nian plus a small perturbation (the 'environment'). The 
square of the scalar product of the initial and final states 
defines the quantum fidelity of the system (Loschmidt 
echo) and has been the object of intense study in re- 
cent years. Jalabert and Pastawski [2] have proven that, 
for perturbations that are classically weak but quantum- 
mechanically strong, the fidelity decay rate only depends 
on the classical Lyapunov exponent of the unperturbed 
system. This universal behavior was later corroborated 
by numerical simulations 0i0 and experiments For 
weaker perturbations, the decay rate is still exponential, 
but perturbation-dependent (Fermi golden rule regime). 
For still weaker perturbations, the decay is Gaussian 

An equivalent approach to the Loschmidt echo was pro- 
posed earlier by Peres 0. In order to study the sepa- 
ration of classical trajectories, it is customary to com- 
pare the evolution of two slightly different initial con- 
ditions. Peres noted that one could just as well com- 
pare the same initial condition evolving in two slightly 
different Hamiltonians, an unperturbed one Hq and a 
perturbed one H = Hq + 5H. The fidelity at time t 
is then defined as the square of the scalar product of 
the wave functions evolving with Hq and H respectively: 



F(t) = KV'.ffoWl ^Hit))\ 2 ■ The latter approach is the one 
adopted throughout the present paper. 

Virtually all theoretical investigations of the 
Loschmidt echo consider one-particle systems evolving 
in a given (usually chaotic) Hamiltonian. The aim of 
the present work is to explore the more realistic case 
of a system of many interacting particles, particularly 
electrons. In order to obtain a tractable model, we shall 
assume that the electrons interact via the electrostatic 
mean field, their dynamics being described by a set 
of one-dimensional (ID) hydrodynamic equations. As 
many experimental studies on quantum information 
involve the manipulation of charged particles, our 
approach may shed some light on the robustness of such 
systems against perturbations from the environment. 

Model. — The physical properties of our model are best 
illustrated by considering its classical counterpart, the so- 
called 'cold plasma' model 0,0|- In the latter, the elec- 
tron population is described by a phase-space distribu- 
tion function of the type: fix, v, t) = n(x, t) 5[v— u(x, t)], 
where n and u are respectively the electron density and 
average velocity, and 6 denotes the Dirac delta function. 
The support of such a distribution function in the 2D 
phase space (x,v) is a ID curve defined by the relation 
v = u(x,t). The electron distribution evolves accord- 
ing to the (collisionless) Vlasov equation. The ions are 
motionless with uniform equilibrium density no and pe- 
riodic boundary conditions (with —L/2 < x < L/2) are 
assumed for all variables. 

As long as there is no particle overtaking in the phase 
space, each position x corresponds to a well-defined ve- 
locity u(x, t). In this case, the Vlasov equation can be re- 
duced to a closed set of pressureless hydrodynamic equa- 
tions: 
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where — e an m are respectively the electron charge and 
mass, and 4>(x, t) is the electric potential obeying Pois- 
son's equation 
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When particles overtake each other, the function x — > 
u(x) becomes multivalued. The above hydrodynamic 
description then breaks down, although the microscopic 
Vlasov model remains valid. 

A trivial stationary solution of Eqs. is given 

by n — no, u = 0. The electron dynamics can be ex- 
cited by modulating the initial velocity: u(x, t — + ) = 
Vocos(kox), with kg = 2n/L, which is equivalent to 
applying an instantaneous electric field at time t = 0. 
The classical dynamics is determined by a single dimcn- 
sionless parameter, namely the normalized wave num- 
ber of the initial perturbation Kq = koVo/tOp, where 
ojp = (e 2 no/m£o) 1 / 2 is the electron plasma frequency. 
Note that Kq can be viewed as the ratio of kinetic to 
potential (electric) energy. It can be shown that two 
different regimes exist. When Kq < 1, electric repul- 
sion dominates, so that the electrons never overtake each 
other in the phase space. In this case, the hydrodynamic 
equations Q"@ can be solved analytically and the so- 
lution displays nonlinear oscillations at the plasma fre- 
quency. When Kq > 1, the analytical solution breaks 
down and the dynamics must be described by the micro- 
scopic Vlasov equation. In the extreme case Kq ^> 1, the 
dynamics becomes again integrable, because it reduces 
to that of free-streaming electrons. For moderate val- 
ues of Kq (but still larger than unity), the electrons can 
be alternately free streaming and trapped by the self- 
consistent (SC) potential. This regime corresponds to 
the formation of complex vortices in the phase space and 
leads to a chaotic dynamics, as was pointed out for the 
similar scenario of nonlinear Landau damping |IJ . In this 
work, we will be mainly interested in the chaotic regime 
and use the value Kq = 2. 

Quantum corrections to the hydrodynamic equations 
0-© were previously derived [l(| ■ For fermions at zero 
temperature (a case that is relevant to electrons in met- 
als), the momentum conservation equation should be 
modified as follows: 



dt dx m d x 2m 2 dx V x /n J mn d x 



(4) 



The second term on the right-hand side is the Bohm po- 
tential: this is a dispersive term that prevents the break- 
down of the quantum hydrodynamics even when Kq > 1 
The third term is the Fermi pressure, which in ID 
can be written as: Pf/Pq = (n/no) 3 , where the equilib- 
rium pressure is given by the usual formula, Po = ^uqEf- 
Ep is the Fermi energy computed with the equilibrium 
density. 

The continuity equation JQ) and the quantum momen- 
tum equation @ can be written in the form of a sin- 
gle nonlinear Schrodinger equation by introducing the ef- 
fective wave function ^(x, t) = y/n(x, t) exp (iS(x,t)/K), 
where S(x, t) is defined according to the relation mu — 
d x S, and n = \^\ 2 . The wave function <3> obeys the equa- 
tion 
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Equation (JSJ) with Poisson's equation Q constitute the 
mathematical model used throughout this Letter. The 
equilibrium Hamiltonian Hq is time dependent, as it de- 
pends self-consistently on the wave function, but con- 
serves both the total mass and the total energy. The ini- 
tial condition is analog to the classical one described in 
the preceding paragraphs and can be easily derived from 
the velocity perturbation u{x, t = + ) by using the rela- 
tion between S and u. Two more dimensionless parame- 
ters (in addition to Kq) intervene in the quantum model: 
(i) the normalized Planck constant h — haj p /mV 2 , which 
measures the importance of quantum effects; and (ii) the 
normalized Fermi velocity vf/Vq. The latter affects very 
little the results (provided it is not too large), and will 
be fixed at vf/Vq = 0.1 in the forthcoming simulations. 
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FIG. 1: Frequency spectrum of the potential energy, for an 
unperturbed evolution with Ko — 2 and h — 0.05. 

The numerical solution of Eq. © is obtained through a 
splitting scheme that separates the kinetic and potential 
parts of the Hamiltonian. Derivatives are computed with 
centered differences. The resulting algorithm is second- 
order accurate both in space and time. 

Results. — First, we characterize the spectral properties 
of the unperturbed Hamiltonian. We consider a case with 
Kq = 2 and h = 0.05 and plot in Fig. 1 the frequency 
spectrum of the time history of the electrostatic energy. 
The spectrum is broad and virtually flat in the range 
< uj < iuj p . The dynamics is therefore sufficiently 
irregular to enable us to compare our results to those 
obtained for a single-particle chaotic Hamiltonian. 

In order to study the behavior of the quantum fi- 
delity, we need to compare the evolution of $ ob- 
tained with the unperturbed and perturbed Hamiltoni- 
ans. We use a static perturbation consisting of a sum 
of a large number of uncorrelated waves: SH(x)/mV 2 = 
e X)j=Af X cos(kjX + otj), where e is the amplitude of the 
perturbation, kj — j (2ir/L), and the ay are random 
phases. The wave number spectrum of the perturbation 
(i.e. the values of iV m ; n and ./V max ) affect only very weakly 
the behavior of the fidelity: therefore, we will focus our 
analysis on the dependence of the fidelity on the ampli- 
tude e. 

A typical result for the quantum fidelity is presented 
in Fig. 2. Contrarily to most single-particle cases, the 
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FIG. 2: Fidelity decay for Kq = 2, h = 0.05, and perturba- FIG. 3: Evolution of the symmetry E(i) for the same case as 



tion e = 10 



in Fie. 2. 



fidelity does not decay exponentially. Instead, it remains 
equal to unity until a critical time r c , after which it de- 
cays abruptly within a few units of u^ 1 . This behavior 
is generic and was observed for all set of parameters that 
were studied, provided the dynamics is sufficiently irreg- 
ular. Numerical tests showed that the value of the critical 
time is independent on the time step. However, as the 
underlying dynamics is chaotic, evolutions computed us- 
ing different time steps will inevitably diverge for long 
times. Therefore, the details of the evolution for t ^S> t c 
are not quantitatively meaningful, and simply indicate 
that the fidelity has dropped to very small values. 

The observed drop in the quantum fidelity is related 
to a sudden symmetry breaking of the wave function. 
Indeed, the evolution equations ©-© for the unper- 
turbed Hamiltonian are invariant under the transforma- 
tion x — —x. If the initial condition is an even func- 
tion of x, this symmetry is thus preserved in time, i.e. 
ty(x,t) = x, t), Vt. But the perturbation SH pos- 
sesses no particular symmetry, and one would expect that 
the symmetry of the initial state quickly deteriorates. 
The symmetry properties can be conveniently measured 
by the following quantity: 
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which is equal to unity when ^f(x,t) = \f r (— a:, £). The 
evolution of S(t) for the perturbed Hamiltonian is plotted 
in Fig. 3. The drop of the quantum fidelity happens 
virtually at the same time as the breaking of the wave 
function symmetry. This behavior is also generic across 
a wide range of parameters. 

We further investigated the dependence of the criti- 
cal time t c on the perturbation amplitude e, for vari- 
ous values of the normalized Planck constant h. The 
critical time is defined as the time at which the fi- 
delity has dropped to 10% of its maximum value, i.e. 
F(t c ) — 0.1. Figure 4 shows that, for small values of 
h, t c depends logarithmically on the perturbation am- 
plitude, i.e. r c ~ — to In e, with ui p t Q ~ 4.3 (this is the 
straight line depicted in Fig. 4). This logarithmic de- 
pendence appears to be universal (both in slope and ab- 



solute value), at least for the values of Kq and vf/Vq 
adopted in these runs. For larger values of Planck's con- 
stant (h > 0.2), this behavior is less neat, particularly 
for small perturbations. 
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FIG. 4: Critical time r c (in units of lu~ ) versus perturbation 
amplitude e, for h = 0.05 (stars), h = 0.025 (diamonds), and 
h = 0.0125 (triangles). The solid line represents the curve 
t c ~ —to In e, with u p to = 4.3. 



A similar pattern was observed for a chaotic quantum 
map m that case, the fidelity stays equal to unity 
until a critical time, after which it starts to decay expo- 
nentially at the classical Lyapunov rate. No sudden drop 
was observed, as is the case for our simulations. 

In order to better evaluate the impact of the SC field, 
we performed some simulations where all nonlinear terms 
have been suppressed in Eqs. The first nonlin- 

earity comes from the Fermi pressure and can be removed 
simply by setting Ep = 0. The SC nonlinearity comes 
from the fact that the electric potential depends on the 
wave function through the electron density n = |^| 2 . To 
remove this nonlinearity, we define the electron density 
independently of 'F, as the sum of traveling plane waves: 

where 



n = n cxt = n [i + a 2_^ j=1 cos(fc J x - u p t + a'j)], 
S is the amplitude of the density fluctuations, and the 
a'j are random phases. This definition is plugged into 
Poisson's equation to yield the electric potential. As the 
resonances of the plane waves overlap in phase space, the 
resulting (time-depend ent) Hamiltonian Hq is likely to 
display chaotic regions [l^] ■ The fidelity decay is studied 
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by perturbing the Hamiltonian in the same way as in the 
SC case. 

In Fig. 5 we plot the quantum fidelity for 5 = 0.5, 
h = 0.025, and several values of the perturbation e. The 
fidelity decay is exponential and begins at t = 0. The 
decay rate is approximately proportional to the square of 
the perturbation, which shows that we are in the so-called 
Fermi golden rule regime |3l 1 1 11]- However, contrarily to 
Ref. [llj, no plateau was observed for short times. 
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FIG. 5: Fidelity decay for a purely external Hamiltonian: e = 
10" 3 (stars); e = 2xl0~ 3 (triangles); e = 5xl(T 4 (diamonds). 
Time is rescaled to the square of the perturbation e. 
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FIG. 6: Fidelity decay for a mixed external and SC Hamil- 
tonian, for (3 = (solid line); /3 = 0.01 (dotted); /3 = 0.03 
(dashed); /3 = 0.1 (dot-dashed); = 0.3 (dot-dot-dashed). 

Finally, we studied a case where both the SC and the 
external fields are present. This is accomplished by defin- 
ing the electron density as: n(x, t) — n cxt + /3(| V E'| 2 — no). 
By varying (3 and 5, we can move continuously from a 
purely 'external' regime (/3 = 0) to a purely SC one 
(S = 0, (3 = 1). We concentrate on a case with 5 = 0.5, 
h = 0.025, and e = 10~ 3 , and vary the value of f3 (Fig. 6). 
For short times, the fidelity decays exponentially with the 
same rate as in the purely 'external' case. Subsequently, 
the decay becomes faster, and for large values of (3 we 
almost recover the abrupt drop of Fig. 2. 

Discussion. — The present work is a first attempt at 
studying quantum fidelity decay in a system of electrons 
interacting through their SC electric field. Our numer- 
ical results show that the quantum fidelity can display 
a rapid decrease. Such effect is probably related to the 
fact that the unperturbed Hamiltonian Hq depends on 
the wave function. When the perturbation SH induces a 
small change in 'J, Hq is itself modified, which in turns 
affects "J, and so on. Because of such nonlinear loop, the 
perturbed and unperturbed solutions can diverge very 
fast (typically, within a few ujp 1 ). In contrast, for the 
single-particle dynamics, Hq is fixed and the solutions 
only diverge because of the small perturbation SH. 

In summary, it appears that the 'natural' response of 
a SC system to a perturbation is a sudden drop of the 
fidelity after a quiescent period rather than an exponen- 
tial decay. Further studies will be needed to understand 
whether these results extend to more complex models, 
going beyond the mean-field approach adopted here. 

We acknowledge fruitful discussions with R. A. Jal- 
abert. 



[1] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

[2] R. A. Jalabert and H. M. Pastawski, Phys. Rev. Lett. 

86, 2490 (2001). 
[3] Ph. Jacquod, P. G. Silvestrov, and C. W. J. Beenakker, 

Phys. Rev. E 64, 055203(R) (2001) 
[4] F. M. Cucchietti et al., Phys. Rev. E 65, 046209(2002). 
[5] H. M. Pastawski et. al, Physica A 283, 166 (2000); M. F. 

Andersen, A. Kaplan, and N. Davidson, Phys. Rev. Lett. 

90, 023001 (2003). 
[6] A. Peres, Phys. Rev. A 30, 1610 (1984). 
[7] J. Dawson, Phys. Rev. 113, 383 (1959); G. Kalman, Ann. 



Phys. 10, 1 (1960). 
[8] P. Bertrand et al., J. Plasma Phys. 23, 401 (1980). 
[9] F. Valentini et al., Phys. Rev. E 71, 017402 (2005). 
[10] G. Manfredi and F. Haas, Phys. Rev. B 64, 075316 

(2001) . 

[11] G. Benenti and G. Casati, Phys. Rev. E 65, 066205 

(2002) . 

[12] A. Macor, F. Doveil, and Y. Elskens, Phys. Rev. Lett. 
95, 264102 (2005). 



